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ABSTRACT 

We explore the rich globular cluster (GC) system of the nearby Sa galaxy M104, the 
"Sombrero" (NGC 4594), using archive WFPC2 data. The GC colour distribution is 
found to be bimodal at the > 99% confidence level, with peaks at (V— J) = 0.96±0.03 
and 1.21 ± 0.03. The inferred metallicities are very similar to those of globular clusters 
in our Galaxy and M31. However the Sombrero reveals a much enhanced number of 
red (metal-rich) GCs compared to other well-studied spirals. Because the Sombrero 
is dominated by a huge bulge and only has a modest disk, we associate the two 
sub-populations with the halo and bulge components respectively. Thus our analysis 
supports the view that the metal-rich GCs in spirals are associated with the bulge 
rather than with the disk. The Sombrero GCs have typical effective (half-light) radii 
of ^2 pc with the red ones being ~30% smaller than blue ones. We identify many 
similarites between the Sombrero's GC system and those of both late type spirals and 
early-type galaxies. Thus both the GC system and the Hubble type of the Sombrero 
galaxy appear to be intermediate in their nature. 
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1 INTRODUCTION 

Observations with the Hubble Space Telescope have, in many 
ways, revolutionized our understanding of globular cluster 
systems (GCSs) in external galaxies. However, ellipticals 
and, to some extent, SO galaxies have received by far the 
largest amount of attention, primarily because early-type 
galaxies tend to have much richer GCSs and suffer less from 
internal extinction problems than spirals. As a result, our 
knowledge of the GCSs in spiral galaxies is still limited to a 
handful of galaxies, including the Milky Way and M31. 

Perhaps the most conspicuous difference between the 
GCSs of ellipticals and spirals is the much higher specific fre- 
quencies (Sn = number of globular clusters per unit galaxy 
luminosity, Harris & van den Bergh 1981) of ellipticals. Spi- 
ral galaxies typically have Sn Ja 1, while nor mal elliptical s 
have Sn ~ 2 — 5 and cDs often reach Sn ~ 15 (Harris 1991). 
Part of this discrepancy may be removed if the Sn of spi- 
rals is normalized only to the spheroi dal component ra ther 
than to the total galaxy luminosity (Cote et al. 2000J), al- 



* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc. under NASA contract No. NAS5-26555. 



though this approach becomes problematic for galaxies like 
the LMC and M33 which lack an obvious spheroidal com- 
ponent but do possess significant population s of globular 
clusters (e.g. Forbes et al. 2000) . McLaughlin (119991) has ar- 



gued that the high specific frequencies of giant ellipticals 
may be explained if the number of globular clusters is in- 
stead normalised to the total galaxy mass, including hot 
X-ray emitting gas. In this way he finds a constant globular 
cluster formation efficiency with respect to the total initial 
gas mass available, but the remaining problem is then to 
explain why the field star formation efficency in the giant 
ellipticals is lower than that of GCs. 

In addition to the Local Group galaxies, some of the 
best studied spirals are M81 and the two edge-on spirals 
NGC 4565 and NGC 5907. Although M81 is nearby, studies 
of its GCs are complicated by internal extinction in M81 and 
by its large angular size which causes confusion problems 
with foreground stars. However, available data indicate that 
the G CS of M81 is very similar to t hat of M31 and the M ilky 
Way (]Perelmuter & Racine 1995[ |Schroder et al. 200l|). In 



NGC 4565 and NGC 5907, internal extinction problems are 
minimized because of the edge-on orientation of these two 
galaxies. Again, the GCSs of these two galaxies (type Sb 
and Sc, respectively) appear to be quite similar to those of 
the Milky Way and M31 in terms of richness and spatial 
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structure , with specific frequencie s of Sn 
al. 199^; [Kissler-Patig et al. 1999J ). 



0.5 (Fleming et 



It has been known for many years that the Milky Way 
globulars can be divided into (at least) two populat ions 



based on t heir metall i cities and spa tial distributions (Kin- 
man 19E9J; Zinn 1985; Minniti 1995), where the metal-rich 



population appears to be associated with either the thick 
disk or the bulge component. Similar conclusio ns have been 



reac hed for the GCS of the Andromeda galaxy ( seitzer et al 
1996 ; J iblonka et al 



1998 



Barmby et al. 2000). In contrast 



Kissler-Patig et al. fll999j ) suggested that, since NGC 4565 
and NGC 5907 contain roughly equal numbers of globular 
clusters but have very different bulge to thin-disk and thick- 
disk to thin-disk ratios, these components of spiral galaxies 
may not have a significant influence on the building of glob- 
ular cluster systems around such galaxies. However, metal- 
poor and metal-rich clusters were not addressed separately 
in the study of NGC 4565 and NGC 5907 and indeed the 
total numbers of GCs detected in these galaxies (40 and 25 
respectively) were too small to permit such an analysis. 

Within the last few years it has been firmly established 
that the GCSs of many, if not most, early-t ype galaxies 
also ex hibit bimodal metallicity distributions (|Gebhardt & 
Kissler- 



Patig 199S; 



Kundu & Whitmore 2001; larsen et al 



2001). The spatial distribution of red (metal-rich) clusters 
tends to be more centrally concentrated than that of the 
blue (metal-poor) ones. One of the more surprising results 
of recent HST imaging has been the discovery that the av- 
erage sizes of metal-poor and metal-rich GCs in early-type 
galaxies differ by ~ 30 %, with the metal-poor clu s ters being 
systemat ically larger ( Kundu fc Whitmore 199£ ; |Kundu et 
al. 1994 |P uzia et aL 1999| ; [Larsen fc Brodie 200C ). The sizes 
of globular clusters in the Milky Way also show a correla- 
tion with galactocentric distance, in the sens e that clusters 



far f rom the centre are systematically larger ( van den Bergh 
1994 J). j This appears to be in contrast to the situation in el- 
liptical galaxies, at least in NGC 4472 and M87 where n o 
such trend is seen (Kundu et al. 1999; Puzia et al. 1999). 



Furthermore, the size difference is appar ently present over a 
large range of galactocentric distances (Larsen et al. 2001), 



so at least in some ellipticals it seems difficult to explain 
it as simply due to differences in the radial distributions of 
metal-rich and metal-poor globular clusters. 

The "Sombrero" galaxy (M104 = NGC 4594) repre- 
sents an important intermediate case between early-type 
galaxies and the Sb/Sc type spirals mentioned above. Al- 
though it clearly contains a star forming disk and a conspic- 
uous dust lane, its integrated luminosity is dominated by a 
huge bulge/halo component. It is one of the near est Sa-type 
gala xies , located at a distance of only 8.7 Mpc (Ford et al 



1996n7jThe first study of globul ar clu sters in the Sombrero 



was carried out by Wakamatsu (1977) who identified a very 
rich GCS with an estimated total of about 2000 GCs. Har- 



ris et al. (1984) identified 1200 ± 100 globular clusters in the 



Sombrero, corresponding to a specific frequency of 3± 1 with 
respect to the bulge/halo luminosity alone. A sim ilar num- 
ber of GCs was found by Bridges & Hanes (1992) who also 
estimated a mean metallicity of [Fe/H] = —0.81 ±0.05 based 
on B — V colours. This mean m etallic ity was confirmed spec- 



troscopically by Bridges et al. (1997) who obtained a mean 
[Fe/H] = -0.70 ±0.3 for a sample of 34 GCs. G roun d-based 
CCD imaging of the GCS by Forbes et al. (|l997|) found 




Figure 1. The three WFPC2 pointings on the Sombrero galaxy 
used for this study, superposed on an image from the Digitized 
Sky Survey. 



some tentative evidence for two populations of GCs. How- 
ever these results were uncertain due to contamination from 
foreground stars and background galaxies in their ground- 
based images. Tw o pop ulations were confirmed by Gebhardt 
& Kissler-Patig (1999) in their study of 50 galaxies using 
archive WFPC2 data. However they did not comment on 
the Sombrero galaxy explicitly. 

Here we present a new investigation of the globular clus- 
ter system of the Sombrero galaxy, using HST archive images 
of three fields. 



2 DATA REDUCTION 

Basic information about the datasets is listed in Table Py and 
the pointings are illustrated in Fig. pL The central pointing 
has the PC chip centred on the nucleus of the Sombrero 
galaxy and consists of 3 x 400 s exposures in the F547M 
band and 3 x 350 s in F814W. The WFPC2 field reaches out 
to a distance of about 4 kpc from the nucleus. The two halo 
pointings were both observed with the F606W and F814W 
filters and are located at 4'6 (11 kpc projected) and 8'3 
(21 kpc projected) from the centre of the Sombrero. Un- 
fortunately, only one exposure of 800 s in each filter was 
available for the inner halo pointing, which makes cosmic 
ray removal a more difficult task. For the outer pointing we 
use 2 exposures of 2100 s each in each filter, far deeper than 
the two other pointings. In addition, we use F606W and 
F814W data for a reference field at 227' from the Sombrero 
to estimate the number of foreground/background sources. 
For the pointings where more than one exposure was 
available in each filter, the individual exposures in each band 
were combined into a single image using the IRAFP] task IM- 

T IRAF is distributed by the National Optical Astronomical Ob- 
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Table 1. WFPC2 archive images of the Sombrero galaxy used in this paper. 



Field 


PID 


PI 


Filters 


Exposures 


Centre 


5512 


Faber 


F547M + F814W 


3 x 400 + 3 x 350 s 


Inner Halo 


5091 


Groth 


F606W + F814W 


800 + 800 s 


Outer Halo 


5369 


Griffiths 


F606W + F814W 


2 x 2100 + 2 x 2100 s 


Reference 


7909 


Casertano 


F606W + F814W 


500 + 400, 500 s 



COMBINE. For the central pointing, background subtraction 
was then performed by first removing point sources and then 
median filtering the resulting images as described in Larsen 
& Brodie ( 200C ) . The median-filtered images were then sub- 
tracted from the original images. 

The region around the dust lane and the high surface 
brightness areas near the nucleus were masked out before 
further analysis. Input object lists for photometry were then 
produced by running the daofind task in the daophot 
package within IRAF on both the V and I band images and 
matching the two resulting coordinate lists. As an additional 
selection criterion, the background noise was measured di- 
rectly on the images in an annulus around each object and 
only objects with a S/N > 3 in both V and I within an 
aperture radius of 2 pixels were accepted. In this way we ac- 
counted for the varying background noise, especially within 
the central pointing. For the inner halo pointing where the 
individual exposures suffered from quite substantial num- 
bers of cosmic ray (CR) events, the object lists were in- 
spected manually to make sure that all sources for which 
photometry was subsequently obtained would be real ob- 
jects. The requirement that objects were detected both in V 
and I turned out to be quite efficient in eliminating CR hits 
from the object lists. However, a few GC candidates were 
manually removed from the object lists because of CR hits 
that would have affected the photometry. 

The V— I colours were measured through an r — 2 pix- 
els apeljfcUXe using thp phot taglr in nAPiPWOT while a larger 
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r — 3 plixels aperture was used for V magnitudes to reduce Clayton fc Mathis 1989J) 



Figure 2. Colour-magnitude diagram for objects in the WFPC2 
chips on the central pointing. The dashed lines mark the bound- 
ary of the region within which globular cluster candidates were 
selected. Typical V — I errors are indicated by the error bars. 



sizes are twice as large. For V — I the aperture corrections 
are 0.030 and 0.120 mag for the WF and PC, changing by 
no more than ~ 0.01 mag for any reasonable cluster sizes. 
Finally, a correcti on for Galactic fore ground extinction of 
Av = 0.17 mag ( Schlegel et al. 1998J ) was applie d to the 
photometry, correspond ing to E(V — I) = 0.086 (Cardelli 



systematic errors on the aperture corrections resulting from 
finite object sizes. The use of a smaller aperture for colours 
is justified because the V and I aperture c orrections nearly 



can cel out when colour indice s are formed (Holtzman et al 



1996 ; L arsen fc Brodie 200C ) . The photometry was trans- 
formed to standard V, I magnitudes using the procedure 
in Holtzman et al. (1995). Globular clusters are expected 
to be resolved on HST images at the distance of M104 so 
we determined aperture corrections to the Holtzman et al. 
r = 0'.'5 reference ap erture by convolving King profiles with 
the TinyTim PSF ( JKrist fc Hook 1997J ) and carrying out 
aperture photometry on the resulting images. Because of the 
different image scales on the WF and PC chips, the aperture 
corrections will also be different. We found aperture correc- 
tions in V of —0.15 mag and —0.55 mag for the WF and 
PC, respectively. These corrections assume clusters with ef- 
fective (half-light) radii R c s ~ 3 pc. For larger clusters our 
magnitudes will be systematically too faint, by ~ 0.2 mag 
on the WF and ~ 0.3 mag on the PC chip if the cluster 



servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc. under contract with the National 
Science Foundation 



In addition to the photometry, we also measured sizes 
for individ ual objects in the HST images using the ishape 
algorithm (Larsen 1999a). ISHAPE convolves the PSF (in this 
case generated by TinyTim) with a series of different model 
profiles, adjusting the size (FWHM) of the model until the 
best possible match to the observed profile is obtained. In 
the case of HST images, is hape also convolves with the 
WFPC2 "diffusion kernel" (JKrist fc Hook 1997| ). For the 
model profiles we adopted King ( P-962Q profiles with a con- 
centration parameter (ratio of tidal vs. core radius) of 30. 
Fitting the concentration parameter and linear size simulta- 
neously would require better signal-to-noise and/or angular 
resolution than what is currently availab le for the Sombrero. 
However, as discussed in Larsen (1999a), the effective radii 



derived from the ishape fits are not very sensitive to a par- 
ticular choice of model profile as long as the sources have 
intrinsic sizes smaller than or comparable to that of the PSF. 



3 RESULTS 

3.1 The central pointing 

We first consider the data for the central pointing, as this is 
where the largest number of globular clusters was found. 
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Figure 3. Completeness functions for the WF frames (solid line) 
and PC frame (dashed line), determined from artificial object 
experiments. 



Fig. shows the reddening corrected colour-magnitude 
diagram for all objects detected in the central pointing 
down to V = 25. The globular cluster sequence can be 
clearly seen, extending from V ~ 20 down to V ~ 24 
and with V — I colours between 0.8 and 1.4. Objects out- 
side this colour range are expected to be mainly background 
sources. We selected globular cluster (GC) candidates as ob- 
jects within the region marked by the dashed lines in Fig. H 
i.e. 0.75 < (V-I)o < 1.45 and V < 24. This results in the 
identification of 151 GC candidates, listed in Table 0. If we 
convert the colour interval to a metallicity range using the 
calibration of Kissler-Patig et al. ( 1998 ), then it corresponds 
to -2.05 < [Fe/H] < +0.24, but we note that these limits 
should only be taken as approximate since they are beyond 
the actual calibrated range. 

Completeness tests for the photometry were carried out 
by adding artificial objects to the images and redoing the 
photometry to see how many of the artificial objects were 
recovered. Because the completeness functions will depend 
on object size, we generated the artificial objects by con- 
volving the TinyTim PSF with King profiles with effective 
(half-light) radii of 3 pc. Since the GC candidates do not 
show a strong concentration towards the centre of the galaxy 
and we are masking out the regions near the dustlane, the 
artificial objects were simply distributed at random within 
each frame. The recovery fractions are shown as a function 
of magnitude in Fig. H for the WF (solid line) and PC frames 
(dashed line). From these tests we estimate the average 50% 
completeness level to be at V ~ 24.5 in the WF frames and 
at V ~ 23.5 in the PC frame. The brighter 50% limit in 
the PC frame is mainly due to the higher background level 
there. We thus reach well below the expected turn-over of the 
globular clus ter luminosity function (GCLF) at My ~ —7.5 
(Harris 1991), corresponding to V ~ 22.4. 



The number of contaminating objects can be estimated 
from the comparison field. This field is located so far from 
the Sombrero galaxy (227' ~ 575 kpc) that no globular 
clusters are expected here. We find about 8 objects in this 
field resembling globular clusters by their colours and mor- 
phological appearance. Some of the objects in the compari- 
son field are unresolved by ishape (presumably foreground 



Figure 4. Luminosity function for globular cluster candidates 
in the NGC 4594 central pointing. The solid line shows the best 
fitting t 5 function, M™ = 22.10 and a t = 0.81. Dashed lines: t 5 
functions with the turn-over shifted by ±0.15 mag (the estimated 
1<t error). Dotted-dashed line: at = 0.91. 



stars) while others have apparent sizes comparable to globu- 
lar clusters. The latter are most likely background galaxies. 



3.1.1 Colour distribution 

A glance at Fig. ti suggests two peaks in the (V — I)g colour 
distr ibution. Applying a KMM test (Ashman, Bird & Zepi 



1994) to the sample of GC candidates, we find that the 



colour distribution is bimodal at the 99.9% confidence level 
with peaks at (V — I)o — 0.96 and 1.21. These V — I colours 
correspond to metallicities of [Fe/H] = — 1.4 and [Fe/Hl = 
—0.54 when using the relation of Kissler-Patig et al. (1998). 
The KMM test itself does not provide error estimates for 
the peak colours but from experiments with different selec- 
tion criteria for GC candidates we find that the peak colours 
are probably accurate to about ±0.03 mag, corresponding 
to ±0.1 in [Fe/H]. The KMM test assigns 45% of the GC 
candidates in the central pointing to the blue (metal-poor) 
peak and 55% to the red one. 



3.1.2 Luminosity function 

The luminosity function (LF) for GC candidates in the cen- 
tral pointing (corrected for Galactic foreground extinction) 
is shown in Fig. M. The LF clearly has a peak at V ~ 22, 
although the histogram in Fig. H has more faint objects than 
would be expected for a "standard" Gaussian LF. The field 
around the Sombrero galaxy has many background galax- 
ies and some of the faintest objects may be contaminants, 
but because the halo fields (see below) and the compari- 
son field do not contain the same large numbers of faint 
objects within the GC colour range, part of the departure 
f rom a Gaussian LF may well be real. As noted e.g. by Seeker 
(1992), there are other analytical functions which actually 



match the faint end of the GCLF in the Milky Way and M31 
galaxies significantly better than a Gaussian, notably the 'ts' 
function. The £5 and Gaussian functions are hardly distin- 
guishable above the turn-over, but for magnitudes fainter 
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than the turnover a £g function has larger numbers of ob- 
jects than a Gaussian. 

We carried out £5 function fits to the data in Fig. U us- 
ing a maximum-likelihood algorithm. A correction for com- 
pleteness was applied, using the completeness function de- 
rived from the artificial object tests, and contaminating ob- 
jects were statistically subtracted from the GC sample us- 
ing the comparison field. Objects in the PC frame were ex- 
cluded from the fits because of the brighter completeness 
limit there. Fitting a £5 function to all data in the WF chips 
down to V = 24, we found a peak at m v ° = 22.28 ± 0.12 
and a dispersion of at — 1.02 ± 0.12. Note that the disper- 
sions of th ( at) and Gaus sian (a g ) functions are related as 



10 



1.29 at (Seeker 1992). Applying a Kolmogorov-Smirnov 



(K-S) test (e.g. Lindgren 1962), the above fit is only ac- 
cepted at the 82% confidence level, indicating that the £5 
function does not provide an incredibly good fit to the ob- 
served GCLF in the Sombrero over the full magnitude range. 
Without the completeness correction the peak would be at 
m v ° = 22.22 ± 0.12, while we get m v ° = 22.46 ± 0.14 if 
both the completeness and contamination corrections are 
omitted. Thus, the correction for incompleteness only has a 
small effect on the turn-over magnitude, while the contami- 
nation correction has a more pronounced effect. 

If the fit is restricted to a narrower magnitude range 
then the effects of contamination, incompleteness and de- 
viations from the assumed £5 shape of the GCLF may be 
further reduced. If we adopt a magnitude limit of V = 23.5 
instead of 24.0 then the resulting £5 fit has a turn-over at 
m v ° = 22.05 ±0.10 and a t = 0.81 ±0.10, corresponding to a 
Gaussian a g = 1.05. The K-S test now accepts the fit at the 
99.8% confidence level, indicating no statistically significant 
deviations from a £5 function for this brighter subsample. 
Leaving out the completeness correction has virtually no ef- 
fect, changing m v ° by only 0.01 mag to 22.04 ± 0.10. This 
is not unexpected, considering that the applied magnitude 
limit is ~ 1 mag brighter than the 50% completeness limit. 
Again, the contamination correction has a stronger effect 
and without it the turn-over would be at m v = 22.13±0.12. 
For an even brighter magnitude cut-off one expects the fits 
to become more uncertain as the sample size decreases and 
the dispersion becomes less well constrained, especially if 
the dispersion and turn-over are fitted simultaneously. Nev- 
ertheless, we find that the turn-over magnitude and dis- 
persion of the £5 function remain stable for a magnitude 
cut-off as bright as 22.0, even though the formal errors be- 
come larger. For cut-offs at V — 23.0, 22.5 and 22.0 we find 
m v ° = 22.10 ± 0.12, 22.10 ± 0.18 and 22.08 ± 0.24, respec- 
tively. 

Another possibility is to keep the dispersion fixed and 
only allow the turn-over to vary. If we adopt a fixed a t — 1.0, 
corres pondi ng to a Gaussian ac — 1.3 (e.g. Kundu & Whit- 
more (2001)) and perform a fit down to V = 23.5 then we 
get a somewhat fainter turn-over at m™ = 22.15 ± 0.13. 
However, the K-S test only accepts this fit at the 66% confi- 
dence level, indicating a much worse match to the data than 
for the two-parameter fit. 

We thus adopt m™ = 22.10±0.15 as our best estimate 
of the reddening-corrected V-band turn-over magnitude for 
the GCLF of the Sombrero, where the ±0.15 uncertainty is 
based on the la errors returned by the maximum-likelihood 
fits as well as the scatter in the various fits. This corre- 



1 


1 1 1 '1 






- 


± 


1 — — 1 




- 


++ + 




- 


± 
± 

± 
+ 


± + 

+ + + + + 

± + ± 

+ -* t + + 

+ + + + + + + + + + + 


+ 
4 


+ _ 



0.6 



0.8 



2.0 



Figure 5. Cluster half-light radii measured on the F547M im- 
ages (central pointing) as a function of (V — I)o colour. The error 
bar indicates the typical ran dom e rror on individual size measure- 
ments, as estimated in Sect 3.1.3. 



sponds to an absolute value of My = —7.60 ± 0. 15 fo r 
a distance modulus of 29.7. For comparison, Seeker (1992) 
found My° = -7.29 ±0.13 and My° = -7.51 ±0.15 for the 
Milky Way and M31. Note that the errors on these numbers 
are quite similar to our formal errors, consistent with the fact 
that a similar number of clust ers a re fitted. More recently, 
Barmby, Huchra and Brodie (2001) found a V-band turn- 
over ofm£° = 16.84 for M31, or My° = -7.6. Ferrarese et 
al. (200C) point out that a weighted mean to GCLF turn- 



over magnitudes for galaxies in the Virgo and Fornax gives 
My = —7.60 with a systematic difference of 0.5 m ag be - 
tween the two galaxy clusters. Kundu & Whitmore (2001) 
favour a value of My = —7.41 based on a sample of 28 el- 
liptical galaxies, but do not include any of the Fornax galax- 
ies. Thus, our turn-over of —7.60 for the Sombrero may be 
somewhat on the bright side, but it is not clear how much of 
the scatter in reported turn-over magnitudes from various 
literature sources is intrinsic, and how much of it is due to 
different measurement techniques etc. The dispersion mea- 
sured for the Sombrero GCLF (a t = 0.81L0.10 or a g = 1.05) 
is formally somewhat narrower than the value in the Milky 
Way, but agrees well with the one in M31 (a g — 1.42 and 
1.06, respectively, Seeker 1992). 

Overall, we conclude that the GCLF of the Sombrero 
seems quite similar to that in most other well-studied galax- 
ies. However, it is worth noting that the relatively small 
number of clusters does not provide very strong constraints 
on the exact form of the GCLF, a situation which can only 
be remedied by obtaining a larger sample of GCs. 



3.1.3 Globular cluster sizes 

Fig. U shows cluster sizes as a function of (V — 7)o colour. 
The sizes were measured by ishape on V frames, excluding 
objects fainter than V — 23.5 in order to avoid too low S/N 
objects where the size information would be uncertain. 

As can be seen from the figure, the data for GCs in 
the Sombrero show a correlation between cluster size and 
(V — I)o colour similar to that in ellipticals and SO galax- 
ies (Kundu fc Whitmore 1998; Puzia ct al. 199£ ; Larsen et 
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0.75<(V-I) <1.05, V<23.5 




Figure 6. Distributions of half-light radii for blue (top) and red 
(bottom) clusters in M104. On the average, blue clusters are larger 
than red ones. Hatched and outlined histograms are for sizes mea- 
sured on F547M and F814W band images, respectively. 




Figure 7. Comparison of object sizes for objects brighter than 
V = 23.5, measured on F547M and M814W images. The solid 
line is a fit to the data points while the dashed line represents a 
1:1 relation. 



al. 2001). The size difference between blue and red clusters 
(dividing at (V — I)o — 1.05) is seen somewhat more clearly 
in Fig. M where the size distributions for blue (top) and red 
(bottom) clusters are illustrated. In Fig. pi the hatched and 
outlined histograms represent sizes measured on F547M and 
F814W band images, respectively. The sizes measured on the 
F814W images are generally somewhat larger (see below), 
but the size difference between red and blue GCs is clearly 
seen in both cases. The mean effective radii (defined as the 
radius within which half of the total cluster luminosity is 
contained) of red and blue clusters are 7? e ff = 1.61 pc and 
Rett = 2.09 pc, measured on the F547M images. 

We can also compare with cluster sizes measured on the 
PC frame. These should be more accurate, given the better 
resolution of the PC chip. For the 9 red and 9 blue cluster 
candidates on the PC chip with size information, the average 
sizes are 1.48 and 2.54 pc. Of course, small number statistics 
play a significant role here and the clusters near the nucleus 



might also have different physical sizes from those further 
out. In any case, the size difference between blue and red 
clusters is confirmed by both measurements. 

The accuracy of the measured cluster sizes can be fur- 
ther checked by comparing size measurements on the F547M 
and F814W images. Fig. M compares sizes measured in the 
two bands for objects brighter than V = 23.5, the same mag- 
nitude limit adopted for the cluster size measurements. The 
dashed line represents a 1:1 relation between the two sets of 
size measurements while the solid line is a least-squares fit 
to the data. As can be seen from the figure, the sizes mea- 
sured on the F814W exposures are on the average ~ 0.75 pc 
larger than those measured on the F547M images. It is worth 
noting, though, that one WF camera pixel corresponds to 
a linear scale of 4.2 pc at the adopted distance of the Som- 
brero. The 0.75 pc offset could be due to small changes in 
the telescope focus ("breathing"), minor offsets between the 
combined exposures or inaccuracies in the PSF modelling 
by TinyTim. Until these effects are better understood, the 
absolute values of the sizes in Fig. B and Fig. M should not 
be taken too literally, but size differences are much more ro- 
bust. Apart from the 0.75 pc offset, the two sets of size mea- 
surements agree quite well. The standard deviation around 
the fit indicated by the solid line in Fig. m is 0.7 pc, indicat- 
ing that individual cluster sizes are accurate at about this 
level, not counting possible systematic effects. 



3.2 The halo pointings 

The colour-magnitude diagrams for the two halo fields are 
shown in Fig. H Obviously, the number of GC candidates in 
the halo fields are much smaller than in the central pointing, 
but these fields nevertheless serve as a very useful compari- 
son. Using the same selection criteria as for the central field, 
the inner and outer halo fields contain 31 and 11 GC can- 
didates, respectively (Table H). While most of the objects 
in the inner halo field are probably true globular clusters, 
we expect a significant fraction of the objects in the outer 
halo field to be contaminants, bearing in mind that about 8 
foreground/background objects are expected within the GC 
colour range. 

The colour-magnitude diagrams in Fig. H are not quite 
as strikingly bimodal as the corresponding plot for the cen- 
tral field (Fig. 0), and bimodality is only detected at the 
40% confidence level by a KMM test when combining data 
for the two halo pointings. Nevertheless, the peaks found by 
the KMM test are at (V-J)o = 0.97 and (V-J)o = 1.19, 
nearly identical to those in the central pointing. 60% of the 
objects in the halo pointings are assigned to the blue peak 
by the KMM test, compared to 45% in the central point- 
ing, providing a hint that the ratio of blue (metal-poor) to 
red (metal-rich) clusters increases somewhat as a function 
of galactocentric distance. 



3.3 Radial trends 

Data for all the three pointings have been combined to pro- 
duce Fig. |9| which shows cluster sizes as a function of pro- 
jected galactocentric distance R s . Red ((V— I)o > 1.05) and 
blue ((V — 7)o < 1.05) clusters are shown with different sym- 
bols. We have not attempted to remove contaminants from 
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Figure 8. Colour-magnitude diagram for objects in the WFPC2 
chips in the two halo pointings. 



Fig. H. Although some of these will be foreground stars and 
thus appear point-like, others may be background objects 
with finite sizes that could be in the apparent GC size range. 
In addition to the contamination problem, observed radial 
trends in the cluster sizes may also be affected by the fact 
that the different exposures are taken through different pass- 
bands, with different exposure ti mes a nd at different epochs. 
In particular, Holtzman et al. (L999) noted that the HST 
point spread function can exhibit significant temporal vari- 
ation due to the so-called "breathing" of the HST secondary 
mirror. 

With these precautions in mind, the plot shows little, 



if any, correlation between cluster size and R g , at least for 
R s ^ 15 kpc where most of the objects may still be expected 
to be globular clusters. Indeed, the average effective radii of 
objects in the three pointings are 1.77, 1.99 and 2.02 pc for 
the central, inner and outer halo fields respectively. 

The three HST pointings also allow us to investigate 
the radial density profile of the GCS as shown in Fig. nol 
The central pointing has here been divided into three bins, 
plotting data for the PC chip separately and dividing the 
WF chip data at 3 kpc or about l!2. The two halo point- 
ings each contribute with one data point. The total surface 
density of blue+red clusters is shown with error bars corre- 
sponding to Poisson statistics. We also plot data for red and 
blue clusters separately, although the error bars are omitted 
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Figure 9. GC sizes as a function of projected galactocentric 
distance in kpc. Plus (+) markers indicate red ((V — I)o > 1.05) 
clusters, diamonds (o) blue clusters. 
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Figure 10. The surface density profile of the Sombrero GCS. 
The two outer points are for each of the halo pointings. The two 
lines represent analytical fits to an i? 1 ' 4 law and a power law as 
discussed in the text. One arcmin equals 2.5 kpc. 



for clarity. The data in Fig. hcj have been corrected for fore- 
ground/background objects by subtracting 4 red and 4 blue 
objects from each pointing. 

The azimuthal coverage of our data is clearly limited, 
but our four outermost data points for the surface density of 
globular clusters n evert heless agree quite well with results by 
Bridges & Ha nes (|1992 ) and with the GC counts quoted by 
Harris et al. (1984h. The earlier ground-based studies were 



unable to probe the central few arcseconds of the Sombrero 
GCS, covered by the PC camera. It should, however, be 
pointed out that the PC data may still be uncertain because 
of possible obscuration from the dust lane which covers ~ 
50% of the PC chip, even though we have attempted to mask 
out the most prominent parts of the dust lane. 

We fitted the GC surface density ctgc as a function of 
galactocentric radius R g with two different model profiles: 
A single power-law fit of the form 



net 
CTGC OC R g 



(1) 
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yields an exponent of a = —1.25 ± 0.20, somewhat steeper 
than the value —0.75 ± 0.15 found by Bridges and Hanes 



(1992) for the central 4' but consider ably s hallower than the 
— 1.82 value found by Harris et al. (1984) for the outer re- 



gions. The actual profile is likely to be better represente d as 
a composite of two power-laws ( Bridges fc Hanes 1992] ), al- 
though our relatively sparse coverage does not allow a more 
detailed investigation of this issue. However, the data do ap- 
pear to be better fitted by a de Vaucouleurs R ' profile of 
the form 



log ctgc = olR\ /a + P 



(2) 



with a = -2.00 ± 0.18 and /3 = 3.39 ± 0.23. This is only 
slightl y sha llower than the value a — —2.12 found by Harris 
et al. ( 1984 )- 

The total number of GCs in the Sombrero can be esti- 
mated by integration of (pi). Adopting an outer radius of 20', 
we thus obtain a total of 1150 clusters which is nearly a fac- 
tor of two smaller than t he ~ 2000 GCs usually thought to 



1992 



exist in the Sombrero (Harris et al. 1984; Bridges & Hanes 



lorbes, Grillmair & Smith 1997). However, the un- 



certainties on the fitted a and /3 values translate to an un- 
certainty of about ±50% on the total number of clusters. 
The number of clusters also depends on the completeness 
correction, as well as on the adopted outer radius of the sys- 
tem. If this is taken to be 30' instead of 20' then the total 
number of GCs is ~ 1290. We have not corrected for incom- 
plete sampling of the GC population due to the magnitude 
cut-off at V=24. Integrating a £5 function with m™ = 22.1 
and at — 0.81 over the V = 19.0 — 24.0 range, one finds 
that only about 4% of the globular clusters would fall out- 
side this range. Alternatively, one can compare directly with 
the observed Milky Way GCLF: For a distance modulus of 
29.7 our cut-off at V - 24 corresponds t o My = -5 .7. Of 
the 144 GCs in the McMaster catalogue ( Harris 1996 ) with 
tabulated My values, 25 GCs or 17% fall below this limit. 
Thus, we estimate that we may have lost between 4% and 
17% of the GCs in the Sombrero because of the magnitude 
limit. 

The RC3 catalogue lists a total reddening-corrected B 
magnitude of 8.38 and (B — V)o = 0.84 for the Sombrero, 
leading to an absolute V band magnitude of My = —22.16 
for distance modulus 29.7. 1150 clusters then correspond to 
a specific frequency of Sn = 1.6±0.8 (where the ±0.8 comes 
from the ±50% estimated above), in reasonable agreement 
with earlier results but somewhat on the low side. Estimates 
of the bulge vs. total luminosity for the Sombr ero range 
between 0.7 3 ( [Baggett, Baggett fc Anderson 1998| ) and 0.85 
(Kent 198£). If we adopt 0.8 as a compromise, the bulge has 



My = —21.92 and the specific frequency with respect to the 
bulge alone becomes S'jv(bulge) = 2.0 ± 1.0. If as many as 
17% of the GCs are fainter than V = 24 that would increase 
SW (total) to around 1.9 ± 0.9 and SW (bulge) to 2.3 ± 1.2. 

If the red and blue cluster populations are fitted sepa- 
rately we obtain a R = -2.07 ±0.19 and a B - -1.91 ±0.21, 
again indicating that the spatial distribution of blue clusters 
is probably more extended than that of the red ones. The 
a values can be converted to effective radii (containing half 
the total number of clusters) for the de Vaucouleurs profiles, 
yielding i? cff jfl = 6.'7±2!6 and R cfi >B = 9.'2±3.'9 for the red 
and blue GCs respectively, assuming that the density pro- 
file continues to infinite radius. If we instead adopt an outer 



radius of 20' as above, which is probably a more reason- 
able assumption, then the respective effective radii become 
R cBM = 4.'5 ± l.'O (11.4 kpc)and R eSiB = 5'A ± l.'l (13.7 
kpc) for the red and blue GCs. 



4 DISCUSSION 

We have clearly detected two sub-populati ons of GCs in 
the Sombrero galaxy, hinted at in Forbes et al. ( 1997 ). Thus, 
the Sombrero is like most large galaxies in this regard. The 
blue sub-population, with (V — 7)o = 0.96 is quite simi- 
lar in colour to the halo GCs in our Galaxy and M31, i.e. 
(V-I)o ~ 0.92 ( |Barmby et al. 2000[ ). It is also remarkably 
similar to the mean value for a sam ple of early-type galaxies 
studied by Forbes & Forte fl200l| ), i.e. (V - 7) = 0.954 ± 
0.008. The metallicity of [Fe/H] = —1.4 corresponding to the 
blue peak is also similar to that of metal-poor halo clusters 
i n our Galaxy ([Fe /H] = -1.59) and M31 ([Fe/H] = -1.40) 
(Forbes et al. 2000). It is therefore likely that the blue GCs 
in all these systems are old and metal-poor. 

The red GCs have (V - 7) = 1.21 or [Fe/H] = -0.54. 
This is very similar to the metallicities of the "metal-rich" 
sub-population in both our G alaxy and M31, i.e . [Fe/H] = 
—0.55 and —0.58 respectively ( pbrbes et al. 2000 ). So again, 
by analogy with known spirals, we associate the red GCs in 
the Sombrero with a bulge/disk population. 

The relation between host galaxy luminosity and the 
colours of globular cluster subpopulations has been stud- 
ied in a number of recent papers. Forbes et al. (1997) 



found a correlation between host galaxy luminosity and the 
colour/metallicity of the red GCs, but no significant trend 
wit h lum inosity for the blue GCs. More recently, Larsen et 
al. ( 2001 ) found correlations at the 2 — 3<r level with host 
galaxy luminosity for the colours of both GC populations. 
Kundu & Whitmore ([2001 ) and Burgarella, Kissler-Patig & 



Buat (2001) found at best weak correlations between host 
galaxy luminosity and GC colours, while Forbes & Forte 
( [2001 ) detected a correlation between the colour of the red 
GCs and host galaxy central velocity dispersion, but no cor- 
relation for the blue GCs. The Sombrero does not represent 
an outlier with respect to any of the above datasets and 
was, in fact, included in the samples studied by both Larse n 
et al. ( [200 1| ) and Burgarella, Kissler-Patig & Buat ( |200l| ). 
Note, however, that the mean colours of the two GC popula- 
tions q uoted here are slightly different from those in Larsen 
et al. J200l| ) (0.939 and 1.184 instead of the 0.96 and 1.21 
values in this paper), even though the two studies use the 
same dataset. This difference reflects the slightly different 
colour / magnitude cuts adopted in the two analyses, but 
is within our estimated ±0.03 mag error margins. In fact, 
with the GC colours in the present pape r the Sombrero fits 
even better onto the Larsen et al. (2001) relation for galaxy 
luminosity vs. GC colour. 

In spiral galaxies the blue GCs are naturally associated 
with a halo population. In our Galaxy, the metal-poor halo 
GCs are distributed nearly spherically around the galaxy 
centre while the me tal-rich popula t ion has a s omewhat flat- 
tened distribution ( Kinman 195S ; p^inn 198C ). The spatial 
distributions of red and blue clusters in the central Sombrero 
pointing are shown in Fig. O. From this figure it seems that 
the spatial distributions of red and blue GCs are not strik- 



© 0000 RAS, MNRAS 000, 000-000 



Globular clusters in the Sombrero 9 



V-l>1.05 




V-l<1.05 




Figure 11. Spatial distributions of red (top) and blue (bottom) 
clusters. The dust lane and the brightest parts of the disk were 
masked out before photometry, so the absence of objects near the 
disk is not real. The apparent peculiar alignment of blue objects in 
the PC frame is probably just a consequence of the small number 
of detections in this chip. 



tre, implying that the ratio of blue to red clusters increases 
outward. 

An interesting difference between the Sombrero GC sys- 
tem and that of other spiral galaxies is the relatively low ra- 
tio of metal-poor (blue) to metal-rich (red) GCs, i.e. about 
0.8:1 in the central pointing and perhaps slightly higher for 
the GCS as a whole. Using metallicities and galactocentric 
radii fo r GCs in th e Milky Way from the McMaster cat- 
alogue (Harris 199C ) and dividing between metal-rich and 
metal-poor clusters at [Fe/H] = —1.0, the corresponding ra- 
tio within the central 4 kpc of our own Galaxy is ~ 1.2 : 1 . 
In M31 it is ~ 1.0 : 1, using data from Barmby et al. ( 2000 ). 



Further out in the halo of spiral galaxies metal-poor clusters 
dominate, leading to global blue-to-red ratios of about 2:1 
and 3:1 in the Milky Way and M31 respectively. For early 
type galaxies the ratio is almost always >1:1 although red 
GCs also here tend to be more centrally con centrated than 
blue ones (Forbes, Brodie & Grillmair 1997). 



Gi ven that th e stellar bulge-to-disk ratio for Sombrero 
is 6:1 ( Kent 1988J ) and that the overall ratio of metal-poor 
to metal-rich globular clusters is much lower than in spiral 
galaxies with less conspicuous bulges, it seems very unlikely 
that the red GCs are associated with the disk but must 
rather be associated with the bulge. If the red GCs were 
associated with the disk we might expect a number roughly 
similar to that in the Milky Way, i.e. about 50 or so. In 
Section B.S we estimated a total number of around 1150 GCs 



in the Sombrero so if there are roughly equal numbers of red 
and blue GCs then that would imply a total of ~ 600 red 
GCs. Thus the situation in the Sombrero lends support to 
the idea that the metal-rich GCs also in our own and other 
spiral galaxies are more closely associated with the bulge 



than the disk (Harris 1976; Frenk fc White 



798§ 



Minniti 



199J ). 

Another similarity between the Sombrero and other spi- 
rals, as well as early-type galaxies, is the relationship be- 
tween GC colour and physical size. Here we find that the 
blue GCs are ~30% larger than the red ones. This is seen 
i n several early typ e galaxies and also in the Milky Way 
( Larsen et al. 2001 ). It is not currently clear what causes 
this difference, which could be either set up at formation 
e.g. because of different conditions in the globular cluster 
progenitor clouds, or due to ongoing dynamical processes. 
As noted above, the blue GC population does appear to be 
more spatially extended than the red, but the current data 
are too sparse to study radial trends in the cluster sizes in 
detail and e.g. check how the GC sizes correlate with galac- 
tocentric distance. It is also possible that the different GC 
populations are on different orbits, in which case dynami- 
cal destruction processes might also affect them in different 
ways. 



ingly different within the central ~ 4 kpc covered by the 
archive images, but the spatial distribution of the red GCs 
is clearly different from the stellar disk of the Sombrero. 
The limited spatial coverage of WFPC2 makes it difficult 
to draw firm conclusions regarding any differences in the 
distributions of red and blue clusters further out. However, 
there are some indications that the number of blue clusters 
declines more slowly as a function of distance from the cen- 



5 CONCLUDING REMARKS 

Using HST imaging of the central few kpc of the nearby 
Sombrero galaxy we have detected ~150 globular cluster 
candidates. We find a bimodal GC colour distribution with 
peaks at (V—I)o — 0.96 and 1.21, correspo nding to [Fe/H] 
= -1.4 and -0.54 flKissler-Patig et al. 1998"| ). The blue GCs 
have a mean size of 2.09 and the reds 1.61 pc. Including ad- 
ditional data for two halo pointings at projected distances 
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of 11 kpc and 21 kpc from the nucleus of the Sombrero, we 
have studied radial trends in cluster sizes and the relative 
numbers of red and blue clusters. We find that the red clus- 
ters appear to be more centrally concentrated in the galaxy 
than the blue ones, although it would be highly desirable to 
confirm this result with better sample statistics. As an Sa 
galaxy, the Sombrero offers a potential link between the GC 
systems of elliptical galaxies and later type spirals. Indeed 
we find several similarities between the GC system of the 
Sombrero and galaxies of all types. 

• The blue GCs have a mean colour (and hence metallicity) 
that is similar to those seen both in large spirals and ellip- 
ticals. 

• The size-colour relationship is qualitatively the same as 
that seen in the Milky Way and several well-observed ellip- 
ticals. 

• The red GCs have an inferred metallicity that is very close 
to the metal-rich populations of the Milky Way and M31 and 
other large galaxies. 

• Without corrections for incomplete sampling of the GCLF, 
the GC specific frequency normalised to the total galaxy lu- 
minosity is SW (total) = 1.6 ± 0.8, or S'iv(bulgc) = 2.0 ± 1.0 
when normalised to the bulge. Corrections for incomplete 
sampling would increase these numbers by between 4% and 
17%. The specific frequency of the Sombrero is thus similar 
to that of field and group ellipticals. 

• The large number of red GCs compared to other spirals, 
combined with the high bulge/disk ratio in the Sombrero 
and the spatial distribution of the red GCs, suggests that 
this GC sub-population is associated with the bulge (rather 
than the disk) component. 

The proximity and rich GC system of the Sombrero galaxy 
clearly warrant further study. Better HST coverage beyond 
the central few kpc would improve sample statistics and 
size/colour trends with galactocentric radius could be bet- 
ter probed. Spectra with 8 — 10 m class telescopes would 
provide direct metallicity and age estimates for individual 
GCs. Kinematics of the two subpopulations could also be 
explored. 
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Table 2. Data for clusters in the Sombrero. V and (V — I)o have 
been corrected for a foreground extinction of Ay = 0.17 and 
Ely — I) = 0.086 mag. The last column is the estimated effec- 
tive (half-light) cluster radius in pc. Note that sizes for objects 
fainter than V = 23.5 (in parentheses) are uncertain. Prefixes 'C, 
'HI' and 'H2' denote objects in the central, outer and inner halo 
pointings, respectively. 
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Table || (continued) 



ID RA(2000.0) DEC(2000.0) 



V (V-I)o 



ReS 



C-001 



C- 
C- 
C- 

c- 
c- 

C- 
C- 

C- 



-002 
-003 
-004 
-005 
-006 
-007 
-008 
-009 
C-010 
C-011 
C-012 
C-013 
C-014 
C-015 
C-016 
C-017 
C-018 
C-019 
-020 
-021 
-022 
-023 
-024 
-025 
-026 
-027 
-028 
C-029 
C-030 
C-031 
C-032 
C-033 
C-034 



C- 
C- 
C- 
C- 
C- 

c- 
c- 

C- 
C- 



C- 
C- 
C- 

c- 

C- 
C- 
G- 



-035 
-036 
-037 
-038 
-039 
-040 
-041 
C-042 
C-043 



C- 
C- 
C- 

c- 

C- 
C- 



-044 
-045 
-046 
-047 
-048 
-049 
C-050 
C-051 
C-052 
C-053 
C-054 
C-055 
-056 
-057 
-058 
-059 
-060 



C- 

c- 

C- 
C- 
C- 



© 0000 



12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
12: 
RAS, 



40:00.31 
39:59.40 
39:58.69 
39:58.65 
39:58.52 
39:58.45 
39:59.58 
39:59.02 
39:58.81 
39:59.76 
39:58.77 
39:59.33 
39:59.92 
39:58.55 
39:59.15 
39:58.92 
39:59.00 
39:59.66 
40:00.12 
39:59.01 
39:59.58 
39:59.85 
39:59.86 
40:00.59 
40:00.59 
40:00.88 
40:00.73 
40:00.63 
40:01.02 
40:00.18 
40:01.30 
40:01.05 
40:00.10 
40:01.00 
40:01.96 
40:01.90 
40:02.06 
40:01.22 
40:00.79 
40:01.02 
40:02.24 
40:02.49 
40:01.01 
40:02.60 
40:01.19 
40:02.03 
40:02.77 
40:01.68 
40:02.34 
40:01.60 
40:01.51 
40:02.37 
40:02.48 
40:02.43 
40:03.35 
40:03.19 
40:03.04 
40:03.04 
40:02.65 
40:03.48 
MNRAS 000. 



11:37:21.4 
11:37:25.2 
11:37:28.2 
11:37:25.8 
11:37:25.9 
11:37:26.5 
11:37:16.7 
11:37:17.8 
11:37:18.7 
11:37:11.8 
11:37:17.6 
11:37:12.0 
11:37:07.7 
11:37:16.4 
11:37:12.3 
11:37:13.6 
11:37:13.0 
11:37:08.4 
11:37:04.6 
11:37:11.1 
11:37:06.9 
11:36:41.3 
11:36:36.0 
11:36:55.4 
11:36:55.4 
11:37:00.6 
11:36:49.6 
11:36:45.3 
11:36:57.5 
11:36:27.4 
11:37:02.4 
11:36:54.1 
11:36:23.0 
11:36:50.8 
11:37:05.5 
11:37:00.0 
11:36:59.3 
11:36:31.0 
11:36:16.4 
11:36:21.9 
11:37:00.8 
11:37:07.3 
11:36:17.3 
11:37:05.9 
11:36:16.3 
11:36:40.2 
11:37:02.0 
11:36:24.9 
11:36:45.6 
11:36:19.5 
11:36:15.5 
11:36:40.5 
11:36:44.0 
11:36:42.3 
11:37:05.2 
11:36:59.4 
11:36:53.9 
11:36:53.9 
11:36:40.9 
11:37:06.0 
000-000 



19.22 
20.88 
21.53 
20.18 
22.85 
20.70 
21.21 
20.14 
22.54 
21.47 
20.11 
20.96 
23.31 
21.52 
22.45 
20.79 
23.08 
23.80 
22.25 
21.96 
23.22 
22.58 
23.92 
21.70 
21.70 
22.60 
21.02 
20.83 
21.94 
22.73 
23.58 
20.38 
22.60 
22.32 
21.70 
22.99 
21.92 
23.87 
22.22 
22.74 
22.55 
21.19 
22.34 
21.37 
21.32 
23.44 
22.65 
22.64 
23.69 
22.51 
21.17 
23.33 
21.46 
22.83 
23.06 
23.83 
23.45 
23.45 
21.14 
22.06 



1.28 
1.15 
1.08 
1.34 
0.90 
1.17 
1.07 
1.09 
1.04 
0.95 
1.03 
0.98 
1.27 
1.17 
0.95 
0.89 
1.26 
1.11 
1.02 
0.78 
1.19 
1.24 
1.39 
1.24 
1.24 
1.31 
0.98 
1.10 
1.01 
1.25 
1.19 
0.97 
1.25 
1.15 
1.24 
1.13 
1.20 
0.83 
1.22 
1.03 
1.35 
1.26 
0.88 
1.00 
0.94 
1.25 
1.11 
1.42 
1.28 
0.90 
1.18 
0.97 
0.91 
1.26 
1.00 
1.34 
1.00 
0.99 
1.17 
1.20 



1.78 
1.55 
1.05 
1.05 
1.92 
1.46 
1.22 
3.06 
3.50 
1.69 
2.54 
1.60 
0.06 
1.72 
3.00 
4.75 
1.34 
(err) 
2.62 
1.22 
0.52 
3.69 

(2.01) 
err 
1.49 
0.84 
1.36 
1.81 
1.62 
1.10 

(2.85) 
1.49 
1.30 
1.36 
0.78 
3.69 
1.62 

(1.17) 
1.56 
2.92 
1.94 
1.10 
5.70 
1.68 
2.14 
5.31 
3.95 

11.66 

(0.78) 
0.78 
1.62 
0.13 
3.37 
0.91 
1.81 

(6.16) 

err 

1.04 

0.97 

4.86 



C-061 
C-062 
C-063 
C-064 
C-065 
C-066 
C-067 
C-068 
C-069 
C-070 
C-071 
C-072 
C-073 
C-074 
C-075 
C-076 
C-077 
C-078 
C-079 
C-080 
C-081 
C-082 
C-083 
C-084 
C-085 
C-086 
C-087 
C-088 
C-089 
C-090 
C-091 
C-092 
C-093 
C-094 
C-095 
C-096 
C-097 
C-098 
C-099 
C-100 
C-101 
C-102 
C-103 
C-104 
C-105 
C-106 
C-107 
C-108 
C-109 
C-110 
C-lll 
C-112 
C-113 
C-114 
C-115 
C-116 
C-117 
C-118 
C-119 
C-120 



12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:40: 
12:39: 
12:40: 
12:39: 
12:39: 



01.74 
02.32 
03.15 
02.58 
03.58 
03.10 
03.30 
03.13 
04.51 
04.81 
03.36 
04.85 
03.56 
04.61 
04.82 
03.39 
05.09 
04.34 
03.96 
05.61 
04.55 
02.34 
04.27 
03.45 
02.98 
02.21 
04.66 
03.32 
02.51 
04.30 
05.30 
06.66 
06.95 
04.85 
04.15 
05.03 
06.80 
05.42 
04.97 
05.82 
05.93 
04.23 
07.33 
04.60 
01.80 
01.42 
01.52 
01.41 
01.83 
00.98 
01.76 
01.81 
01.42 
00.96 
01.66 
00.22 
59.90 
01.12 
59.85 
59.80 



-11:36:09.2 
-11:36:21.6 
-11:36:46.1 
-11:36:27.7 
-11:36:59.3 
-11:36:44.0 
-11:36:29.7 
-11:36:03.1 
-11:36:45.3 
-11:36:52.1 
-11:36:05.9 
-11:36:52.4 
-11:36:09.5 
-11:36:42.1 
-11:36:42.4 
-11:35:56.0 
-11:36:48.7 
-11:37:14.4 
-11:37:19.1 
-11:37:09.6 
-11:37:18.1 
-11:37:49.4 
-11:37:43.2 
-11:37:48.8 
-11:37:53.8 
-11:37:59.4 
-11:37:43.5 
-11:37:59.1 
-11:38:06.8 
-11:38:01.7 
-11:37:57.1 
-11:37:53.0 
-11:37:51.1 
-11:38:05.4 
-11:38:11.1 
-11:38:14.7 
-11:38:03.2 
-11:38:14.3 
-11:38:17.8 
-11:38:12.5 
-11:38:12.9 
-11:38:26.9 
-11:38:06.7 
-11:38:26.7 
-11:38:02.5 
-11:37:50.3 
-11:37:57.6 
-11:37:55.5 
-11:38:20.8 
-11:37:55.2 
-11:38:21.6 
-11:38:26.9 
-11:38:15.3 
-11:38:12.3 
-11:38:35.6 
-11:37:53.0 
-11:37:44.5 
-11:38:30.5 
-11:37:49.0 
-11:37:55.4 



22.18 
21.11 
23.90 
21.10 
23.65 
22.67 
22.19 
20.37 
23.57 
21.35 
23.95 
21.85 
23.56 
21.56 
20.91 
21.36 
21.42 
22.21 
19.88 
22.66 
21.73 
23.10 
23.52 
23.07 
22.75 
19.57 
22.43 
23.31 
23.23 
22.52 
22.16 
21.41 
20.93 
20.21 
21.56 
21.92 
23.39 
23.04 
23.95 
22.43 
22.32 
22.16 
22.69 
21.50 
21.88 
23.85 
22.12 
22.47 
22.30 
20.16 
21.36 
22.11 
22.04 
21.87 
23.48 
19.98 
22.24 
23.14 
22.02 
21.98 



1.08 
0.95 
1.02 
0.94 
1.04 
1.07 
1.22 
1.11 
0.87 
0.92 
0.81 
0.99 
1.22 
1.06 
0.90 
0.90 
1.20 
1.18 
1.12 
1.39 
1.14 
1.28 
1.38 
1.10 
1.10 
1.27 
1.13 
1.33 
1.09 
1.14 
0.96 
0.95 
1.21 
1.05 
1.05 
0.93 
0.81 
0.99 
1.12 
1.22 
1.31 
1.30 
0.84 
1.23 
1.32 
1.15 
1.21 
1.12 
1.26 
1.23 
0.96 
1.10 
1.13 
1.10 
1.03 
0.96 
1.22 
0.95 
0.96 
1.23 



1.17 
0.45 

(0.32) 
1.43 

(1.04) 
1.88 
1.56 
2.20 

(1.68) 

2.27 

(10.89) 

2.72 

(0.84) 
2.20 
2.20 
2.79 
0.91 
0.58 
3.37 
0.78 
1.62 
0.32 

(0.97) 
1.81 
1.17 
2.46 
2.01 
0.39 
1.62 
0.71 
1.49 
1.81 
0.78 
1.30 
0.65 
1.23 
2.79 
3.11 

(0.58) 
0.52 
err 
1.17 
0.91 
0.78 
1.56 

(0.19) 
1.49 
2.85 
2.33 
1.17 
1.68 
1.43 
2.27 
1.49 
12.96 
2.07 
1.10 
3.11 
1.17 
1.49 
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C-121 


12:39:59.68 


-11:37:52.4 


22.20 


1.17 


1.04 


Hl-01 


12:40:20.67 


-11:31:09.2 


23.62 


1.13 


(3.65) 


C-122 


12:40:01.13 


-11:38:44.0 


21.03 


1.15 


2.98 


Hl-02 


12:40:20.31 


-11:31:12.3 


22.32 


1.18 


5.31 


C-123 


12:40:00.74 


-11:38:49.5 


21.49 


0.88 


1.62 


Hl-03 


12:40:20.11 


-11:30:26.5 


21.94 


1.35 


0.65 


C-124 


12:39:58.99 


-11:37:51.2 


23.34 


1.18 


err 


Hl-04 


12:40:21.98 


-11:30:07.4 


22.54 


0.89 


6.03 


C-125 


12:39:59.99 


-11:38:25.7 


21.78 


1.18 


1.62 


Hl-05 


12:40:25.72 


-11:30:60.0 


21.78 


1.00 


2.59 


C-126 


12:40:00.33 


-11:38:37.5 


23.93 


1.18 


(12.96) 


Hl-06 


12:40:23.62 


-11:31:40.4 


23.29 


1.07 


6.35 


C-127 


12:39:59.40 


-11:38:08.0 


21.69 


0.98 


4.34 


Hl-07 


12:40:25.23 


-11:32:03.3 


23.29 


0.86 


14.71 


C-128 


12:39:59.00 


-11:38:00.4 


23.79 


1.13 


(err) 


Hl-08 


12:40:28.03 


-11:31:43.5 


21.25 


0.86 


0.06 


C-129 


12:39:59.88 


-11:38:30.0 


22.53 


1.11 


1.88 


Hl-09 


12:40:23.88 


-11:32:37.2 


23.04 


1.23 


0.45 


C-130 


12:39:59.09 


-11:38:06.5 


22.81 


0.98 


0.71 


Hl-10 


12:40:19.21 


-11:32:07.0 


22.46 


0.97 


2.72 


C-131 


12:40:00.19 


-11:38:51.9 


21.26 


1.26 


1.04 


Hl-11 


12:40:20.67 


-11:32:49.9 


22.21 


0.93 


0.26 


C-132 


12:39:59.21 


-11:38:19.5 


20.91 


1.15 


1.04 


H2-01 


12:40:01.31 


-11:32:17.5 


21.67 


1.13 


2.83 


C-133 


12:39:58.78 


-11:38:09.5 


22.43 


1.35 


0.84 


H2-02 


12:40:02.74 


-11:31:57.5 


23.91 


1.23 


(0.45) 


C-134 


12:39:59.25 


-11:38:25.8 


20.29 


1.24 


1.04 


H2-03 


12:40:02.91 


-11:32:12.8 


23.84 


1.18 


(23.33) 


C-135 


12:39:59.43 


-11:38:35.8 


21.18 


1.22 


1.30 


H2-04 


12:40:03.43 


-11:32:23.1 


23.42 


0.95 


2.98 


C-136 


12:39:57.97 


-11:37:57.7 


21.53 


1.06 


1.49 


H2-05 


12:40:04.24 


-11:32:21.1 


21.71 


1.20 


1.56 


C-137 


12:39:59.33 


-11:38:43.8 


21.10 


1.19 


1.30 


H2-06 


12:40:04.49 


-11:31:59.1 


20.28 


0.92 


0.32 


C-138 


12:39:58.63 


-11:38:21.7 


23.85 


0.84 


(15.55) 


H2-07 


12:40:04.49 


-11:31:57.8 


22.04 


1.41 


0.06 


C-139 


12:39:58.63 


-11:38:21.7 


23.85 


0.84 


(err) 


H2-08 


12:40:04.93 


-11:32:02.2 


22.88 


1.13 


7.45 


C-140 


12:39:58.54 


-11:38:19.6 


21.68 


1.14 


1.94 


H2-09 


12:40:05.55 


-11:32:31.7 


20.45 


0.97 


2.33 


C-141 


12:39:58.94 


-11:38:35.6 


21.99 


1.16 


1.30 


H2-10 


12:40:06.01 


-11:32:15.7 


20.51 


1.01 


0.19 


C-142 


12:39:58.26 


-11:38:18.5 


21.14 


1.08 


1.36 


H2-11 


12:40:06.15 


-11:32:30.2 


23.74 


0.78 


(0.52) 


C-143 


12:39:58.05 


-11:38:17.2 


23.51 


0.96 


(2.92) 


H2-12 


12:40:06.40 


-11:32:38.9 


21.85 


0.82 


0.06 


C-144 


12:39:59.38 


-11:39:01.5 


21.03 


0.92 


2.20 


H2-13 


12:40:06.75 


-11:32:26.9 


21.24 


1.35 


0.71 


C-145 


12:39:58.31 


-11:38:39.5 


20.09 


0.89 


2.40 


H2-14 


12:40:03.55 


-11:33:10.8 


22.76 


0.99 


1.23 


C-146 


12:39:56.98 


-11:38:01.7 


21.68 


0.98 


1.56 


H2-15 


12:40:07.53 


-11:33:10.0 


22.02 


0.80 


0.06 


C-147 


12:39:57.72 


-11:38:28.3 


22.34 


0.84 


4.28 


H2-16 


12:40:05.53 


-11:33:18.0 


23.97 


0.92 


(8.49) 


C-148 


12:39:58.48 


-11:38:58.7 


22.60 


1.08 


1.36 


H2-17 


12:40:06.14 


-11:33:21.3 


22.40 


0.97 


0.84 


C-149 


12:39:58.37 


-11:38:57.6 


21.89 


1.00 


3.11 


H2-18 


12:40:07.24 


-11:33:26.2 


22.49 


1.07 


3.63 


C-150 


12:39:58.16 


-11:38:54.6 


22.11 


0.84 


2.27 


H2-19 


12:40:03.00 


-11:33:30.6 


22.19 


0.97 


1.04 


C-151 


12:39:57.39 


-11:38:32.1 


19.41 


0.90 


16.33 


H2-20 


12:40:06.34 


-11:33:29.9 


21.44 


1.06 


0.97 














H2-21 


12:40:04.39 


-11:33:36.1 


23.75 


0.99 


(17.50) 














H2-22 


12:40:02.07 


-11:33:16.2 


19.13 


1.23 


0.71 














H2-23 


12:40:01.21 


-11:33:12.6 


22.68 


1.03 


9.53 














H2-24 


12:40:01.02 


-11:33:35.9 


23.76 


1.16 


(3.37) 














H2-25 


12:40:00.51 


-11:33:04.9 


23.80 


1.12 


(err) 














H2-26 


12:40:00.45 


-11:34:02.7 


22.87 


1.12 


2.98 














H2-27 


12:40:00.09 


-11:34:05.8 


22.04 


1.14 


3.43 














H2-28 


12:39:59.79 


-11:33:00.5 


23.90 


1.13 


(8.81) 














H2-29 


12:39:58.63 


-11:33:13.7 


22.64 


1.06 


2.07 














H2-30 


12:39:58.04 


-11:33:19.4 


21.14 


0.96 


2.01 














H2-31 


12:39:57.88 


-11:33:55.0 


21.13 


0.84 


0.39 
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